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Abstract 



We investigate the difference between those relativistic models based on 
interpreting a realistic nucleon-nucleon interaction as a perturbation of the 
square of a relativistic mass operator and those models that use the method 
of Kamada and Glockle to construct an equivalent interaction to add to the 
relativistic mass operator. Although both models reproduce the phase shifts 
and binding energy of the corresponding non-relativistic model, they are not 
scattering equivalent. The example of elastic electron-deuteron scattering in 
the one-photon-exchange approximation is used to study the sensitivity of 
three-body observables to these choices. Our conclusion is that the differ- 
ences in the predictions of the two models can be understood in terms of the 
different ways in which the relativistic and non-relativistic S'-matrices are re- 
lated. We argue that the mass squared method is consistent with conventional 
procedures used to fit the Lorentz-invariant cross section as a function of the 
laboratory energy. 
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I. INTRODUCTION 



Realistic relativistic quantum mechanical models of interacting nucleons are important 
for modeling medium and high-energy reactions involving nuclear targets. The simplest 
generalization of the non-relativistic two-body Hamiltonian is 

p , 

h = — + V M = 2 + m2 + V'. (1) 

m 

In this case the mass operator, M, replaces the non-relativistic center of momentum Hamil- 
tonian, h. Bakamjian and Thomas have determined sufficient conditions on V in (|^) so 
M can be interpreted as the mass operator of an equivalent relativistic quantum model. 

In the non-relativistic case, interactions are determined by fitting parameters in the 
interaction so that the solutions of the dynamical equations reproduce the existing bound- 
state and scattering observables. This is a non-trivial task that has led to the construction 
of a number of realistic interactions |H [1]. Schiavilla has repeated this analysis, 
by refitting parameters in an Argonne V14 interaction, to determine a realistic relativistic 
interaction, V. It is tedious to have to refit parameters for every new interaction. 

A direct way to utilize existing non-relativistic interactions, without refitting the inter- 
action, is preferable. Two methods exist that can be used to directly formulate relativistic 
models based on realistic non-relativistic interactions. Historically the first of these meth- 
ods was introduced by Coester, Pieper, and Serduke [H- In their construction an 
interaction is added to the square of the non-interacting mass operator. The interaction 
is chosen so the dynamical equations are equivalent to the non-relativistic equations. The 
result is that the eigenfunctions of the relativistic model can be trivially obtained from those 
of the corresponding non-relativistic model. The second method was introduced by Kamada 
and Glockle P]. This utilizes an S*- matrix preserving unitary transformation, called a scat- 
tering equivalence [|ll|] [Q, that relates the two operators in (|l]). The transformation 
involves a simple rescaling and a renormalization that ensures unitarity. Both methods lead 
to invariant S'-matrices that are identical to the non-relativistic S-matrix. 

In this paper we investigate the difference between the mass-squared method and the 
Kamada-Glockle method by comparing both the two-body dynamics and elastic electron- 
deuteron scattering, which is the simplest three-body reaction. 

The surprising conclusion of our analysis is that even though both methods give the 
same invariant scattering matrix as the corresponding non-relativistic model, they do not 
lead to identical predictions, even at the two-body level. Specifically, these two methods do 
not lead to scattering equivalent theories. 

We also consider differences in the predictions for elastic electron-deuteron form factors 
based on these two models. While normally a change in the two-body interaction would 
lead to a corresponding change in the strong current, in typical calculations the consistency 
of the dynamics and the strong current is normally limited to ensuring current conservation 
and current covariance. 

The comparisons in this paper are limited to understanding the quantitative sensitivity 
of the physical observables to the two different models of the two-body dynamics, keep- 
ing everything else the same. Specifically we compare the calculated observables using 
the deuteron eigenstates constructed from the mass-squared and Kamada-Glockle methods 
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based on the same realistic nucleon-nucleon interaction. In both calculations the current is 
treated in an identical impulse-like approximation, using identical nucleon form factors and 
an identical form of the dynamics. It is impossible to rigorously determine if the difference 
in these three-body observables are differences that would arise if the two methods were 
two-body scattering equivalent, or if they are due to the violations of the scattering equiva- 
lence. However, the result of the above comparison is that there are small but measurable 
differences in the two predictions. The parameter that controls the scale of these differences 
is the same one that controls the scale of the violation of the scattering equivalence. 

For the electron-scattering calculations the coupling to the electron is treated in the one- 
photon-exchange approximation. The relativistic models are formulated in Dirac's point 
form of the dynamics. The hadronic current operators are modeled using the point-from 
spectator approximation [O]. The input nucleon current matrix elements are constructed 



from the Gari-Kriimpelmann parameterization [1^. The specific choice of realistic interac- 
tion, form of dynamics, or nucleon form factors is not central to our analysis. 

The notion of scattering equivalence is necessary to understand the issues involved in 
our analysis. Unitary transformations preserve probabilities in quantum mechanics. The 
conventional interpretation of this observation is that the physical predictions of a quantum 
theory do not depend on the choice of Hilbert space basis. 

There is a more subtle interpretation of this result that is best understood in the context 
of scattering theory. The S-matrix measures the difference between the free and interacting 
Hamiltonian that can be observed asymptotically. Scattering observables are preserved 
under unitary transformations provided the same unitary transformation is used to transform 
both the free and interacting Hamiltonians. 

A class of unitary transformations, called scattering equivalences, have the property that 
if they are applied to the full Hamiltonian without applying them to the free Hamiltonian, 
they give the same scattering matrix elements and bound state observables as the original 
Hamiltonian. They are not meant to transform representations of free particles. In this 
sense they have a different interpretation than a standard unitary transformation. They 
parameterize the insensitivity of the asymptotic S'-matrix observables to certain short-range 
features of the interactions. 

Scattering equivalences are unitary transformations W that satisfy the asymptotic con- 
dition |ni g 



lim ||(/-l^)e-*^°*|<l>)|| = 0. (2) 

t — >dioo 

Examples of scattering equivalences include operators of the form 

(3) 

where C is any compact Hermitian operator. 

In a field like nuclear physics, where model interactions are determined by fitting to S'- 
matrix and bound-state observables, realistic model Hamiltonians can at best be determined 
up to an overall scattering equivalence. The ambiguity in the choice of representation of the 
interaction is particularly confusing in the case of relativistic models, where there are many 
scattering equivalent ways to interpret a given non-relativistic interaction. For example, 
the kernel of a single two-body interaction defines different operators in each of Dirac's fl^ 
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forms of the dynamics. The operators differ by whether they commute with the total hnear 
momentum, the total four velocity, or the light-front components of the four momentum. 
The resulting two-body models are scattering equivalent, but lead to different three-body 
predictions. Even within a single form of dynamics, there are scattering equivalences 



that relate different natural choices of two-body interactions. In all of these cases, the 
different scattering equivalent models have the same two-body bound state and scattering 
observables. They can be made to have the same three-body observables by introducing 
additional three-body interactions |]14|. If, as is customary, three-body interactions are 
not included, then different scattering equivalent two-body models give different three-body 
predictions. 

In what follows a brief description of each model is given. Differences in the 5'-matrix, 
wave functions, and elastic electron-deuteron observables are analyzed. 

II. THE MASS-SQUARED METHOD: 

In this section the mass-squared method is described for the case of two particles with 



equal masses. The case of unequal masses is treated in [17 



Let Vnr be a non-relativistic nucleon-nucleon interaction defined by a kernel of the form 

(p, k, j, cr, /, s I Vnr\p', k\ i\ a\ s') 

= (5(p - p')5jj'5aa'(A;, / , s|K„||A;', s'), (4) 

where k is the magnitude of the relative momentum. The physics is contained in the reduced 
kernel 

{k,l,s\\vlJ\k\l\s'). (5) 
Define a model relativistic mass operator: 

:= 4(fc2 + m^) + 4:mV^2- (6) 
where 1/^2 is a point-form interaction defined by 

(v, /c, j, cr, /, s\Vm'^ I V ', k\ i\ a', s') 

= 6{v -v')6jj>5^^>{k,l,s\\viJk',l',s'). (7) 

The interaction kernel of is identical to the non-relativistic kernel (^. In this expression 
V is the velocity of the non-interacting two-body system, k is related to the invariant mass 
of the non-interacting system by Mq = 4(/c^ -|- m^), and j is the spin of the non- interacting 
two body-system. 

The two-body eigenvalue problem, 

M'\^)=X'\^), (8) 
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is mathematically equivalent to the non-relativistic Schrodinger equation 

HnM) = + Vnrm = -enM) ■ (9) 

The eigenfunctions are identical functions. Because of the relation 

= 4m{Hnr + m), (10) 
the eigenvalues of these equations are related by 

X'^ = {2m-erf = Am{-enr + m). (11) 
It follows that for the deuteron 

6y (^nr 1 



4m 2000 



(12) 



which implies that the relativistic and non-relativistic binding energies are identical to within 
experimental uncertainties. 

The two-body mass operator (^ also leads to the same cross section as in the non- 
relativistic model (P). To see this note that the Lorentz invariant |]18| relativistic cross 
section is given by 

/ da \ 



r 



^|(/+|M-MoK)|V^ 
Vr der 



^^\{f+\M-Mo\t)\'k^^ 
2/c/a;'^'' ' ^' 2k 



(2^)-'i{nK„.-i.>p(|^' 



2k/m^^'' I I /I 2k 



\ \<l I I / I 7 

where the relativistic and non-relativistic velocities are 



^'^''^\{nVnrm^^, (13) 
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and the relativistic and non-relativistic kinematic internal energies are 



^2 

= 2u!, and e^^ = 2m H . (15) 

m 

The calculation in (|l^) shows that the relativistic cross section based on is identical to 
the non-relativistic cross section derived from the non-relativistic Hamiltonian (|^). 

This shows that the relativistic invariant cross section (with the relativistic velocity and 
phase space factor) is identical to the non-relativistic cross section (with the non-relativistic 
velocity and phase space factor). The two-body eigenstates of the non-relativistic Hamilto- 
nian (^) are also eigenstates of the relativistic mass operator (^. 

In this case, because the factors /c^ in equations (H) and (Qj are identified, the scattering 
matrix elements are identified as functions of invariant relative momenta k: 



Srik) = SUk). (16) 

These momenta have different interpretations in the relativistic and non-relativistic cases 
when they are expressed in terms of single nucleon degrees of freedom. Both correspond 
to the same asymptotic observable, the momentum of a single particle in the center of 
momentum frame of the non-interacting system. 

Ignoring the difference in relativistic and non-relativistic binding energy, the mass 
squared model is scattering equivalent to the non-relativistic model. There is no need to 
refit interaction parameters. 



III. THE KAMADA-GLOCKLE INTERACTION 



The Kamada-Glockle method is an alternate method that is used to make relativistic 
quantum models using a realistic nucleon- nucleon model as input. 
The Kamada-Glockle method defines the operator q by 

2Vm2 + p := 2m + — = Fiq). (17) 
m 

Equation (p!7|) can be solved for q{k) and k{q): 

q^ = m(2Vm2 + P - 2m); (18) 

Kamada and Glockle also define h{q) by the equation 

edk = h\q)q'^dq (20) 
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which can be solved to give 



ft^M-- 1 + ^ ■ (21) 



Defining the operator Vkg by 

(v, /c, j, a,l,s\Vkg\v', k', j', a', s') 

= 5(v - v05./5..'^^(^, ^II^LII^', I', ^')^^ (22) 

where (/c, Z, A;', ■s') is the reduced interaction (^), Kamada and Glockle show that 

the operator 



M = 2VWT^ + Vkg - 2m (23) 

has the exact same deuteron binding energy as predicted by the non-relativistic modeL 
In addition, they show that the S-matrix elements predicted in their model and the non- 
relativistic model, ClS db function of 

, o2 

e = 2Vm2 + F - 2m = — , (24) 

m 

are identical. Thus, these two models are scattering equivalent. 

The relativistic wave functions and transition matrix elements can be obtained directly 
from the corresponding non-relativistic quantities: 

cj>^{kJ,s) = ilj^{qik)J,s)/hiq{k)) (25) 

and 

T/(fc, I, s; e; k\ s') := J^fMk^ /, e; q{k')^\ (26) 

This has the advantage that the solution of the relativistic scattering and eigenvalue problems 
do not have to be recalculated. 

With this method, the non-relativistic and relativistic scattering matrix elements are 
identified as identical functions of the non-relativistic and relativistic invariant center of 
momentum energy: 

Sr{e) = Snrie); (27) 

where 

e = q'^/m + 2m = 2\/WTm?. (28) 
Thus, the Kamada- Glockle model is also scattering equivalent to the non-relativistic model. 



7 



IV. DISCUSSION 



Each of the methods discussed above provides a simple prescription for interpreting ex- 
isting reahstic non-relativistic models as equivalent relativistic models. In the mass-squared 
case a local interaction is added to the square of the free mass operator and in the Kamada- 
Glockle case a non-local operator is added to the mass operator. 

Because each model has a relativistic scattering matrix that is equivalent to the non- 
relativistic scattering matrix, one might expect that the two relativistic models would be 
scattering equivalent. This is not true. This is because in the mass-squared case the S- 
matrix elements are identified as functions of k while in the Kamada-Glockle case they are 
identified as function of the invariant energy. 

In order to get a quantitative understanding of the difference in the S-matrix elements 
of these two models, note that both models solve the Schrodinger equation: 

2 J 2 

(- + = ei2|0) (- + V^2M = ei2|0). (29) 

m m 

Although the operators are different, the internal wave functions are the non-relativistic wave 
function as a function of g or respectively. It follows that the corresponding scattering 
matrix elements are related by 

SKcif^ikKa)) = Snrienriq)) = S,n2{krr,2) for k^2 = q. (30) 

The problem is that at point where all of the S-matrix elements agree, we have /c^^ = Q ^ 
kxG- On the other hand, kxG and fc^a represent the same asymptotic observable in the two 
relativistic models, which means the the physics requires 

S^2{k^2) = SKoikKG) for k^2 = kKG (31) 



which contradicts 

To show this explicitly note by definition 

kla = q'a + qV^ra'). (32) 

The point where the mass squared and Kamada-Glockle S'-matrices agree corresponds to 
km2 = q. Substituting /c^2 for in the above equation gives the following relation between 
k^2 and kxG- 

kla = kl2{l + kl2/Am'). (33) 

Since kKG > k^2 at the point where the S-matrix elements agree, this implies that the 
KG cross sections fall off slower than the mass-squared cross sections when considered as 
functions of k. 

There is nothing theoretically wrong with either of these methods; nevertheless they are 
inconsistent with each other and cannot both be consistent with experiment. 



The key to resolving this problem was originally noted by Breit |19]. He correctly pointed 



out that the laboratory cross section is a Lorentz invariant function of a Lorentz invariant 
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variable, the laboratory energy. Specifically, the laboratory energy eiab can be expressed in 
a manifestly invariant manner by: 

Pi ■ P2 = -eiabm = -m^ - (34) 

or 

(^iab = '^ hm, k"^ = —{eiab - m). (35) 

m 2 

This indicates that the experimentally determined S-matrix and cross sections are Lorentz 
invariant functions of the invariant eiab or equivalently k"^. The key point is that there is no 
such thing as the experimentally determined non-relativistic cross section. 

"Non-relativistic" potentials are determined by adjusting the parameters of the potential 
until the scattering eigenstates of 

H=— + V (36) 

m 

reproduce the experimentally measured S^k"^) = S{eiab)- The analysis of the mass-squared 
method in equation (|1^) shows that the non-relativistic and relativistic formulas for the 
cross section lead to the same invariant quantity when k"^ is identified with the asymptotic 
relative momentum in the center of momentum frame. 

If this same data were used to determine the Kamada-Glockle interaction, it would be 
appropriate to determine the interaction Vkg directly by comparing to the Lorentz invariant 
cross section as a function of k"^, as was done by Schiavilla f^. The identification ([T7|) defines 
q by 

g2 = m{2^kl^2+m^ - 2m) (37) 

rather than by k^'^- This would lead to a "non-relativistic" model with a different high q 
behavior. 

In the case of the Kamada-Glockle method, because the measured invariant cross section 
is given as a function of kj^Q rather than g^, the fitting should be done as a function of k"^, 
as in the mass-squared method, with the corresponding q dependence determined by the 
identification 

q^ = m{2^m? + kli^Q — 2m) = m(2'^w? + —{eiab — "m) — 2m). (38) 

This suggests that a Kamada-Glockle interaction constructed from a non-relativistic 
model fit to the S'-matrix as a function of laboratory energy will not be consistent with 
experiment at high energies. In practice, the differences are negligible as long as /c^/4m^ is 
small. For momenta where this quantity is large there are questions about the interpretation 
of non-relativistic models. 

Although the wave functions in the two models are not observable, they are shifted in 
the same manner as the S-matrix. The rescaling and renormalization of the wave function 
is responsible for all of the differences in our computations of the electron scattering ob- 
servables. The result is that the Kamada-Glockle wave functions will be spread out more in 
momentum space than the mass-squared deuteron wave functions. 
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The factor k"^ /Am^ in (|33|) defines the scale of the shift of the S'-matrix or wave function 
at a given value of k that results from the two different interpretations of the k variables. 
This factor also defines the scale at which the two models fail to be scattering equivalent 
and the scale on how much the the factor h{q'^) in (^) differs from unity. This suggests 
that the differences in the wave functions of the two models are not simply due to an 
underlying scattering equivalence. Instead, they are completely determined by the difference 
in the momentum dependence of the 5- matrices. In the case of the electron-scattering 
predictions, using these wave functions in the same point-form spectator approximation 
leads to differences in the elastic observables. While this in not a true test of the three-body 
sensitivity, this comparison is consistent with the way that models of the hadronic current 
are tested. 

Figures 1 and 2 show the momentum space s and d components of the deuteron wave 
functions for the mass-squared and Kamada-Glockle models based on the same Argonne 
V18 potential. These are plotted as functions of the variable /c, which is a function of 
invariant mass of two free particles. In the mass-squared case this is also interpreted as the 
non-relativistic relative momentum, while in the Kamada-Glockle method the corresponding 
"non-relativistic" relative momentum, g, is smaller. As discussed above, this tends to spread 
out the momentum space Kamada-Glockle wave functions relative to the mass-squared wave 
functions when plotted as functions of k. This is seen in the plots of the s and d-staXe wave 
deuteron wave functions in Figures 1 and 2. 

Figure 1 shows the reduced s-state wave function and Figure 2 shows the reduced d-state 
wave function. The solid line in both figures corresponds to the mass-squared model while 
the dashed line corresponds to the Kamada-Glockle model. Both figures show the predicted 
shift of the Kamada-Glockle wave function to the right of the mass-squared wave function. 

The shift is due to the relationship between k and q in (^2]) rather than due to the 
factor h{q), which is a compensating factor that simply ensures that the change of variables 
is unitary. Because this shift is minimal when q^ /Am? <^ 1, the two wave functions look 
similar for momenta out to about 5fm~^ [q^ /Arn? ~ .25). 

We conclude that the difference in these two models is due to the difference in how the 
non-relativistic invariant energy is interpreted relativistically. Specifically with the mass- 
squared method, fc^ is identified as a relativistic relative momentum while in the Kamada- 
Glockle method the relativistic and non-relativistic rest energies are identified. These dif- 
ferences explain how the models fail to be scattering equivalent. Nevertheless, both models 
provide reasonable relativistic interpretations of non-relativistic models. They differ on 
scales where the non-relativistic interpretation is in question. 

These differences have implications for the three-body problem. Because current matrix 
elements can couple high and low momentum components of the wave functions, the possi- 
bility of observable differences between these two approaches at momentum transfers below 
l(GeV)^ exists. This is particularly relevant in the point-form spectator approximation, 
where the momentum transfer seen by a constituent nucleon is larger than the momentum 
seen by the deuteron [jT^| 

The calculations shown in Figures 3, 4, and 5 are the deuteron elastic electron scattering 
observables ^(Q^), B{Q^), and T2o(Q^) respectively in the point-form spectator approxi- 
mation In the point-form spectator approximation matrix elements of the transverse 
current. 
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>(0) := J^(0) - q^^^^. (39) 

are identified with matrix elements of tlie one-body operator, 

>(0) = Jm + ^p"(0), (40) 

in tlie Breit frame. All of the other current matrix elements are determined by current 
conservation and current covariance. The one-body currents are determined using Gari- 



Kriimpelmann single-nucleon form factors. The dotted curves are the the result of 
using the Argonne V18 interaction to construct a Kamada-Glockle interaction and solid 
curves are the result of using the same interaction to construct a mass-squared interaction. 

The results show small, but measurable differences, in A{Q'^) at momentum transfers 
near 1 (Gev)^. The observables B{Q'^) exhibit differences in the two predictions that extend 



below the scales where differences are expected. The data in Figure 3 is from refs. |^0[| ||21 
2| m while the data in Figure 3 is from [||] [|3l || 



The results for the tensor polarization, T2o(Q^), are also similar. Differences in the 
predictions of the two methods again extend below = l(GeV)^. The data are from |27 



We conclude that the difference in these two methods is due to differences in the interpre- 
tation of the relative momentum variable that appears in the non-relativistic Hamiltonian. 
The relevant rescaling is given in (^). The variable k"^ /Am? defines a scale where the dif- 
ferent interpretations of the momenta become significant. This is the scale on which the 
two methods fail to be scattering equivalent and is the scale that defines the shift in the 
momentum space wave functions. The mass-squared method is consistent with parameter- 
izations of the Lorentz invariant cross section as a function of the projectile energy in the 
laboratory frame. For a fixed current model, the three-body calculations are sensitive to the 
differences in the wave functions. The model calculations show that these differences can 
have a measurable impact at much lower momentum transfers. To consistently interpret 
the Kamada-Glockle model relativistically, the non-relativistic phase shifts would have to 
be refit. This only impacts the high momentum structure of the model. 

These differences are shown to have measurable implications in elastic electron-deuteron 
scattering. In the realistic model presented, these differences lead to measurable effects in 
the observables -B(Q^) and T2o(Q^) for momentum transfers as low as a few tenths of a 
(GeV)^, far lower than the scales where differences are expected. 

We conclude that the mass-squared and Kamada-Glockle methods define different dy- 
namical models and that the mass-squared method is the one that is consistent with con- 
ventional methods for determining experimental cross sections. Quantitively, the differences 
in these models are small, but they have dynamical consequences at surprisingly low mo- 
mentum transfers in the three-body system. 

This work was supported by the Department of Energy, Nuclear Physics Division, under 
contract DE-FG02-86ER40286. 
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FIGURES 



FIG. 1. The momentum space deuteron s-wave wave function for a relativistic point- 
form dynamics based on the mass squared model (solid line) and the Kamada-Glockle 
model (dashed line). 



FIG. 2. The momentum space deuteron d-wave wave function for a relativistic point- 
form dynamics based on the mass squared model (solid line) and the Kamada-Glockle 
model (dashed line). 



FIG. 3. A{Q'^) for the mass squared interaction (solid line) and the Kamada-Glockle 
model (dashed line). Both calculation use Argonne V18 as the input interaction, Gari- 
Kriimpelmann |]13| form factors, and the point-form spectator approximation. Data is 
from pg PH PU pi. 



FIG. 4. B{Q'^) for the mass squared interaction (solid line) and the Kamada-Glockle 
model (dashed line). Both calculation use Argonne V18 as the input interaction, Gari- 
Kriimpelmann |]l3l form factors, and the point-form spectator approximation. Data is 
from pi pi pi pg. 



FIG. 5 T2o(Q^) for the mass squared interaction (solid line) and the Kamada-Glockle 
model (dashed line). Both calculation use Argonne V18 as the input interaction, Gari- 
Kriimpelmann [|l3l form factors, and the point-form spectator approximation. Data is 
from p^ p8[ PU pO] 0. 
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